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ABSTRACT 


The determination of the optimal distribution of aimpoints is 
examined for weapons that fire fragmenting projectiles against mobile 
targets. A finite difference approximation which reduces the problem 
to a mathematical programming problem is developed. Computational 


considerations for this nonlinear programming problem are discussed. 
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JL. NATURE OF THE PROBLEM 


Consider a situation where an artillery battery or naval gunfire 
support ship desires to fire at a moving land target. If the original 
position, speed, and direction of the target were known, the guns 
could be aimed so that the projectiles would impact at a calculated 
future location of the target. This future location is fairly easily 
calculated in this simple situation. The calculations become more 
difficult whenever any of the factors of original position, speed and 
direction are unknown. In these cases a probabilistic determination 
of future target location becomes necessary. It is then possible to 
state only that the target is more or less likely, in a probabilistic 
sense, to be at a certain location than at another. | 

This paper will examine the problem of determining the optimal 
distribution of aimpoints for a weapon which fires a projectile against 
a moving target whose speed is known but whose original position and 
the direction of movement are not. 

A model for the location of a target under the above circumstances 
was developed by B. Koopman in 1946 (Ref. 1), and is described herein. 
Although formulated in a Naval setting, Koopman's model equally well 
applies to a moving target that has been detected by a forward observer; 
the exact position of the target being unknown, All that is known is 
that the target is more likely to be at a point O than at any other 
point. The target may not be’at 0, however, but only within a short 


distance of 0, all points the same distance r from 0 being equally 





likely. The probability that the target is in an area dA at a distance 
xr from the target is defined as P(xr)dA. Koopman assumed that the 


situation could be approximated by the circular normal distribution, 


(2) = exp {2/20 , (1) 





where < is the variance in the original target location. 
The speed of the target, u, is assumed to be known but the direction 
of the target movement is unknown, all directions being equally likely. 


After t units of time the situation may be as shown in Figure 1. 


dA winaALy 





Figure 1. Entry of target into dA 


Koopman determined the distribution of moving targets about a point O 


to be 





PG ce > exp {-(x°4u*?)/2s%, 1, (3) (2) 


where 1. is the modified Bessel function of the first kind and zero 


order. The graph of P(r,t) for different t is shown in Figure 2, 
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Observe that the probability spreads outward with time so that the 


target is most likely to be in an expanding ring about 0. 
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Figure 2. The distribution of moving targets about a point O. 





IIT. MATHEMATICAL MODELS OF WHAPONS FIRING 
FRAGMENTING PROJECTILES 


It was necessary to examine the mathematical formulation of a 
weapon that fires fragmenting projectiles in order to determine the 
distribution of aimpoints at a moving target. The basic models for 
such weapons usually include four probability distributions. 

The probability that a target located at coordinates Xi9 Vy is 
killed by a round impacting at coordinates x,y is defined as P. (Xx, -*,7,-1)- 
This probability may be considered as the lethality function. The 
lethality function is the conditional probability that a target at “Sec 
Yas is killed given that a round impacts at x,y. It may take several 
forms of which exponential, linear and "cookie-cutter" are familiar. 
Reference 2 discusses lethality functions. 

The distribution of impact points about the point of aim is defined 

as PL(x,y). This distribution is caused by meterological effects on 
the projectile and ballistic effects inherent in the weapons system. 
Some models have incorporated the distribution of impact points about 
the point of aim into the lethality functieam It is then possible to 
define the probability of killing a target located at coordinates Xi9 Vy 
with a round aimed at coordinates Xo» Vae This probability is defined 
as 19 Cavers een) 

The probability distribution of aimpoints for a target is dependent 


upon the nature of the target, its movement and location, and the type 


of weapons system employed. It is defined as P (x9 ,)- 





The probability distribution of the location of the target is 
defined as Pry (4 294.) 

These distributions have been combined to model weapons systems. 
This paper examines three of these models to illustrate the modeling 
techniques and to provide background for examining the problem of 
determining the distribution of aimpoints of a weapon firing fragmenting 
projectiles against a moving target. Although a different notation is 
used in each of the original source documents, a standard notation has 


been adopted for purposes of presentation in this thesis. 


A. THE GROVES MODEL 

This model was developed by Groves in Ref. 3. I1t is a simple model 
that does not consider a distribution of aimpoint; i.e. the aimpoint 
feetaxed. If ce cay eeoe is the lethality function representing the 


probability that a target locatea at coordinates x is killed by a 


os 
round impacting at coordinates x,y, and P,(x,y) is’ the distributroensor 
the impact points; then the probability of killing a target at Xi 9 Vy 
by a round aimed at coordinates x9 Vo is 


© © 
a) oe J P, (x,-*,y,-y)P,(x,y) dx dy. (3) 

The probability that a target survives one round is then He Ogre oe): 

and the probability of surviving N rounds all aimed at the same aimpoint 

is (1-P,(x,-",.7,-¥)) « The probability that the target is killed by 

any of the N rounds is 1-(1-P,(x,-%,55,-7,)) « This expression can then 

be integrated over the area of the target, T, to obtain the expected 


fractional kill, K 


K = fl (1-(1-P,(x,-, 57 ,-¥,)) ae, dy,- (4) 
if 





The substitution of equation (3) into (4) yields 


© 


K= [fQ-Q- fof PCy vyw)PgGay) & 4") ax, oy. (5) 


Pa © 0] om CO 


Groves further reduces this expression to one that is more suitable for 
hand—computing the expected fractional kill. It is noted that the 
simplicity of this model is a result of the fixed aimpoint as seen by 


comparing the Groves model with the Weiss and the Breaux-Mohler models. 


B, THE BREAUX=-MOHLER MODEL 

This model was published by Breaux and Mohler in Ref. 4. The model 
contains all of those distributions discussed earlier. It is used by 
Breaux and Mohler to compute the expected fraction of a target damaged 
by a salvo of fragmenting projectiles all aimed at the same aimpoint. 
They define the probability that a target located at coordinates 
X,, Y, is damaged by a round impacting at coordinates x,y as PL (x,-*,y,-y)- 
P,(x,y) is the density function describing the distribution of impact 
points x,y about the aimpoint x9 Yee It is assumed that all N rounds 
of a salvo are aimed at the same aimpoint. The probability a target 


survives N rounds is 


© oe) 


O- ff PCeexy ny) (ey) ae ay)", (6) 


and the probability of damage over all impact points is 


Lo @) (oe) : 


1 (ae ff PL Gepeayyny)P_(oeyy)ax ay)”, (7) 


Breaux and Mohler assume that the target is distributed over the target 


area, T, as Pra(Xp 904) 9 and the aimpoint itself is a random variable 





distributed as P(x 19.) The expected fraction of the target damaged 


is determined to be 


co co oe) roe) , N., | 
f = JJ J J. | a-(1- les ‘ Py Gey) Py Gry) dx ey Po yyy) 
P(x yg)& oy, oc ady,- 8) 


This expression is reduced using a binomial expansion to. 


N ; fo) Oo co ’ 
Feria per rls f P, (eg-xy¥,-¥)P, (xy) ac ay | 
J= =CO cD —{oo} —(s2) 


Pyle Fy JEgle I a)%, Ve, O, O° (9) 


Breaux and Mohler further reduce equation (8) using Jacobi polynomials 
to produce an expression that can be used for determining the expected 


ieaetional kill. 


C. THE WEISS MODEL 

This model was developed by Weiss in Ref. 5. It will be discussed 
more thoroughly than the Groves and Breaux—Mohler models because the 
model is used later in the paper. Weiss assumes a target of n men that 
is distributed as P(X V4) The probability that a man is in the 
: ax, dy, and 


small area &, a is P(X Vy 


ff Palxpryeqaay = 2. (20) 


The expected number of targets in a small area is nP,, (x, »¥, dx, dy, 
The probability that a round aimed at the aimpoint x Y, will kill a 
target is I See Weiss assumes that N rounds are fired at 


the target and the ith round has an aimpoint Xo 


ees sare The probability 


10 





that a ee at amon survives all N rounds is 
N 
a(x, 4 -¥q) — (1-P; (4-5 5%4-Vas)) 
and the probability that a target is killed is 
1=a(x,-*_39.-7,)- 


The probability that there is a target at x9 Va 


Pap s¥4)(1-a(x,-% 974-9, ; 


The expected number killed (K) is 


a 
tt 


‘ J, J. Py (oego¥g)(I-alyo¥4-¥Q) Og 


CO co c CO 


\ 


CO co =i 
ce ES J. J Ply 1) a Gy, 14, -V) OH, a | 


The expected fraction surviving ? is =, or 


ioe) ce 


n 


On 
i} 


Loe) 


i Pr(x, »¥,Ja(e,-*% 74-7) cx aes 


pao @) pate o) 


If each round has the same aimpoint, then 


oO co : 


N 
$(x_s7,) = J J 1 Caer, J Ma are ie ayy )) COS i 


and using a binomial expansion 


co 


BiG 5), 


ce 


j=0 


11 


BRL LL Mayra ayer) aye] 


co C N - : . 
N,N 5 
J | PralX, 2¥4,) a (4) P(e, TEV) dx, dy. 


(11) 


(12) 


and Ghaceis Kaseled as 


(13) 


1G J J Pr(Xps¥, Oy “nl : Prlx, s¥4)a(X,-%, 594 -¥,) 449,» 


(14) 


(15) 


(16) 


N o roe) : . ; 
N,N 4 
ECG) FL Ppl Gyre va)” PM. (17) 





Assuming that the aimpoints themselves are distributed as NCAEN eer 


the expected fraction of targets surviving is 


1H 
Il 


J be P(x, 9) OC, 97, )O wv > (18) 


Or 
| 


oe) oe) N ro ro) 
a _1)X (N = Sed 
_ J P (x, s¥,) ms 1) () ic _ P(X 94) Pyle xsl 5) 
° dx, dy, dx dy. (19) 


D. COMPARISON OF BREAUX-—MOHLER AND WEISS MODELS 
It is easily shown that the Breaux-—Mohler and Weiss models give the 
same results. Substituting equation (3) into equation (8) gives the 


expected fraction damaged as 


fs Uh J J 1-(1-P (p74) Pp leprry)Py (Xq sg) 4S 4% QV, (20) 


foe) 
T OO a 


Substituting equation (16) into equation (18) gives the expected fraction 


of targets surviving as 
co © co © N 
= _ J 15 CSE ora) J J ipl Coup Mle es, eo) dx, dy, dx, dy(21) 


It is clear that these models are simply complements of each other, 


III, THR DETERMINATION OF NECESSARY AND 


SUFFICIENT CONDITIONS 


The mathematical models described in chapter II can be used to 
determine necessary and sufficient conditions for determining the 


aimpoint distribution that maximizes the expected fraction of the 


12 





target destroyed or minimizes the expected fraction of survivors, 
depending upon the model used. Reference eee mathematical models 
of hit probabilities and the techniques of maximizing salvo kill 
probabilities. Weiss (Ref, 7) attacks this problem using techniques of 
Svesnikov (Ref. 8) and Morse and Kimball (Ref. 9). His approach to the 
problem of maximizing the salvo kill probability (or minimizing the 
fraction of survivors) required the application of the calculus of 
Variations. This application of the calculus of variations can be 
understood by examining the problem considered by Morse and Kimball in 
Ref. 9. Their treatment is presented herein to facilitate the under- 
standing of Weiss! development of the necessary and sufficient conditions 


contained in Appendix A. 


A, THE MORSE AND KIMBALL PROBLEM 
Morse and Kimball (Ref. 9) examined the problem of determining the 
firing pattern that maximizes the probability of at least one hit on a 
target when N rounds are fired in a single salvo. An approximate 
solution is determined for large patterns. 
Morse and Kimball define the probability that a eae will hit 
the small area element dxdy as f(x,y)dxdy. Then 
eon es . 
f f f(x,y)dxdy = N, (22) 
—0 00 
Pere N is the number of rounds fired. The function f(x,y) is considered 
to be the pattern density function. Extensive changes in the hit 
probability can be obtained by changing the firing pattern. The expected 


number of lethal hits on a target located at coordinates x,y is L 1 Geehi)h 


19 





where L is the lethal area of the target. Morse and Kimball approximate 
the expected number of hits on the target using the Poisson probability 
distribution. The probability of at least one hit is L-exp{t(x,y)}; (23) 
which is the probability of destroying the target. The total probability 
of destroying the target is then 


P= , f (ane OH) ) 6 (yy)axay, (24) 


—_O _O 


where a is the probability density for aiming the pattern, usually the 
normal density. The problem then is to determine the function f(x,y) 
which maximizes P subject to equation (22). The problem is maximize 


pal f (ree MOH) ) 6 (ayy )axay, 


—_OO =O 


Oo © 


Subject to { f £(x,y)dxay = N, (25) 


HOS. Sp = (0). 
It a the method of solution of this problem by Morse and Kimball that 
is interesting. They consider a pertubation of the function f(x,y); that 
is, an increase of f(x,y) by a small amount 6 at the point x4 sVys and a 


decrease of f(x,y) by the same amount 5 at x The constraint remains 


A ee 


satisfied while the objective function changes by an amount 


ELEC, y5) 9 ( KLf(%55To) 6 (a 
(c 17° 4 £ (4 974) -e ORES fo (xps¥5) )é dxdy. 


) are > 0; then if 


Suppose f(x, 574) and £(x517, 


aie : -Lf(x | | 
e 7 a) Ns ? Pe! Xx a9,); 


P can be increased by a 6 > 0. Conversely, if 


lia Ge | —-Lf(x 
e ( ree, 1¥,) —< ( 202) f (xy¥5)s 


14 





then P can be increased by 6 <0. Hence, for that function f(x,y) 


which maximizes P, it must be the case that 
colle cr yc | -Lf(x,,y | 
e bE (xy » Vt (qo¥4) -e (X51 2)f (x59 5) (26) 


for all points x,y where f(x,y) > 0. For all such points 


cu a. a6 = 0 (27) 
Now, instead of x,,y, a point X39 where f(x,y)=0 is considered. If 
£(x4 54) is decreased by 6 (which now must be positive) and f(x; 9¥3) is 
increased by 6, then the increase in the objective function is 

(£(x5975) - oO Tye (ay,)) 5 axay, (28) 


which equals (£5 (x50¥5) - c) 6 dxdy. This implies a positive increase in 
P is possible if £ (x3 53) >c. Hence, f(x,y) cannot equal 0 unless 


£ (xy) Sc, It is seen that the solution is 


§(x,y) 


a (an) 
1 é 


(29) 


= 0 cin £ (xy) G- 


The unknown constant c is determined from the constraint 





N = ee Gyr (30) 


ey) 
The set of conditions of equations (29) form a set of necessary and 
sufficient conditions for determining that pattern that maximizes the 
total probability of destroying the target. 

Weiss' analysis of the problem of minimizing the fraction of targets 
surviving a salvo of N rounds is presented in detail in Appendix A. The. 
original paper (Ref. 7) is greatly condensed and the entire development is 
presented in this paper to facilitate the understanding of the techniques 


msed in vhe solution. 


19 





IV. THE SOLUTION 


A, ANOTHER REPRESENTATION OF THE PROBLEM 


Weiss' coverage problem discussed in Appendix 1 can be written as 


Co © 
Minimize =f f P,(x,,y,) o Ue) as, ay. 3 


Subject to © @ 
JS Palxaey,) o, oy, 


II 
= 


(31) 


Ceo) =o 


where 
[oe @) 


ule) = JS Pax sy.)in aly, y-¥,) ox, 9y,. 


P(x, 334) and P (x, 994) are the target and aimpoint distributions 
described earlier. It is noted that in addition to the constraints 


above, it is trivially true that 
Co 2 
J. J P(x, 9¥4) = 1. (32) 


The solution to this problem is determined in Appendix 1 to be 


Co ® 


=11( x : 
“J. J. Pox s¥4) © si Ve) tna (xy 4 Yay) Oza =C if Pa(x,sy,)?0 (3) 
2G il Eien 


where a(x,- ee, relates to the lethality function. It is not 


Fair t 
possible to simply substitute other than trivial target distributions or 


lethality functions into the necessary and sufficient conditions of 


equations (33) and determine the distribution of aimpoints that 


16 





minimizes the number of survivors. However, a solution to these 
equations may be obtained by using A Meade difference approximation to 
the integrals. This paper assumes that the solution to the finite 
difference approximation converges to the solution of the original 
problem. The proof of this convergence, however, is beyond the scope 
or this paper. It is realized that finite difference techniques may 
guarantee solution for specific values but usually preclude the seeing 
of relationships between equation parameters. 


Equation (31) can be approximated as 


k k 
Minimize < »» Palapye nee), 
x =1 y.=1 
(34) 
k k 
Subject to = oy P(x, »¥,) = N, 
x =1 y,=t 


Pea 
PPR aid O, 
k k 
where u(x, »4) a z : = P (x, »y,)Ina(x,-x, » FeV 5 ) 
Ps 


and the upper limit of summation, k, is the number of increments in the 
approximation. Equations (33) can similiarly be approximated by 


k k 
. a Paltyry4) © te Ntona (xyes 14 Van ee ae 
cat (35) 
So sige P(x, 9¥,,)=0. 


—u(x 


The corresponding finite approximation to equation (32) is 


k k 
A= P(X, V4) = 1. (36) 
x =1 ¥,=1 


It is noted that equations (35) are the derivatives with respect to 


ce sy.) of the objective function of equation (34) when the expression 


17 





for u(x, »¥,) is substituted into the objective function. Equations (45) 


Can be Wwerttemsas 


a6 =Cc fie P(x, »¥,) ee, 


OP, (x_+¥,) = (37) 


ee P(x, 9¥,) = 0 
It is noted that the objective function is convex in coe and 
that equations (35) and (37) are simply the results of Gibbs liam! 
as applied to equations (34), 

A solution to equations (35) can be obtained with a computer by 
dividing the target area into cells sufficiently small for approximation 
of the small areas dx, dy, and. dx dy. To simplify the computer procedure, 


the following notations will be used: 


P. = The expected fraction of the target in cell i; it is 


analogous to the target distribution. 
y. = The proportion of the total rounds aimed at cell j3; it is 


analogous to the aimpoint distribution P(x, 9¥,)3 


ia The probability that a target in cell i is killed by a 
round aimed at cell j; it is analogous to the lethality 
function q(x,-x_»7,-¥,)- 


The problem in equation (34) can then be rewritten as 


k &k 
k oe ye 
Minimize %® = 2 P.e ist jnt J a 
i=! 
k 
Subject to E P, = 1, (38) 
; i=1 ; 
5 1 
Re Be 
jet” 
Pay J) ae =, 0 


loipbs Lemma is discussed in Ref. 10. 
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The solution ‘1s imeterms of the ars and k is the total number of cells. 


Equations (35) and (37) can be rewritten as 


oy; ij e Ale] =C ciate ig (39) 
or, equivalently, 
P.b.. ej ij = ¥ ify. >0, 
ig J (40) 
< ¥ ube Ys =o | 


where Y = -C. 
The problem thus presented is similar to one presented by Professor 
John M. Danskin to his Games of Strategy class at the U. S. Naval 


Postgraduate School in December, 1970. 


B. THE COMPUTER SOLUTION TECHNIQUE 

A solution technique is to divide an area of one Kilometer square 
into cells of twenty meters square. The P.ts are then determined for 
each of the 2,500 cells from the target distribution. To solve the 
problem of firing at targets obeying Koopmans! moving target distribution 
described earlier, values of time of impact and speed of the target are 
inputs and.the P,'s are determined from equation (2) for each cell of 
the simulated target area. ‘The ee are determined from the lethality 
functions for the weapon system being investigated. For an artillery 
weapon it is the case that the Gia = 0 for all cells j more than, 
say, 60 meters from cell i. This means the number of machine calculations 
required for solving the problem in equations (40) would be reduced 


considerably as the sum over j would be limited to the 48 cells 


immediately surrounding cell i and, of course, cell i itself. 


1) 





Tt 2s) netedsthatens Diy = 0 for all j # i; that is, the probability 


that targets outside the aimpoint cell are killed is zero, the solution 


is 
aoe 
y.=o-in— if P.bue iva 
). eee eke 
=e if “Pope. eee 
neal 


yY is determined from the equation 


k 4 P. D.. 
id D se id 
i=l “ij Mf 


which is analyzed in Ref. 10. 
After the values for the Pts have been initialized and the lethality 
function formulas are determined for the cells surrounding each cell i, 
a beginning set of values for each ys is initialized. A possible set of 
y,'s is Velale Voreee= ¥o50079° It is then necessary to determine if that 
initialized set of ye give a minimum value of @ and, if not, in which 
direction to move (i.e. what other vector of y;'s) to produce a lower %. 
An algorithm has been developed by Professor Danskin that will cause 
convergence to a vector of fe that produce a minimum ®, The computer 
program for this algorithm as well as the application of the algorithm 
to the problem discussed in this paper are presented as a separate thesis 


for Professor Danskin by Major Paul T. Zmuida, USA. 


20 





ts 


APPENDIX A 


Weiss desires to minimize 
CO ew 


N 
= J J Pn) i a(x,-X 4 s¥.-¥,,) coe oy, (A-1) 


Subject to 


co 


J J P(X, 174) dx, dy, =a. 
Pas s¥4) 20, 


where $ is the expected fraction of the target surviving, ete hee, is 


the target distribution, and nm q(x 


7 1¥4,-Y,3) is the probability a 


=—_x 
[Rapes il 


target at coordinates x survives a salvo of N rounds when the ith 


4% 


round is aimed at aimpoint Xa Vaile Weiss wishes to minimize ® by a 


peoper choice of x_.,y 


lak Velie 


PeneCy Lune tlOnl To sGellied clchn. taat 


N 
u(x, 74) = -Ln a q(x,-x. 5 Y4-Y a3): 


7s [ala Tyan) oe OC Tay) Ls 


{In Ge MaaNege) Wooo Ine ay a4 Vay | 


I 


N 
a In q Cl Ss aii eeat ine) 
N 
or u(x, 5¥,) gee Ln q (x,-%. 5 9Y47Y aa): (A~2) 


The objective function can then be rewritten as 


(oe) 
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Weiss then makes an approximation that a distribution function for the 

aimpoints can be substituted for the exact knowledge of X49, 80 that 

the number of rounds aimed at area dx dy_ is P,(x ,y_)dx dy... The total 
a*a A*‘ a’‘ae ata 


number of rounds is 


J J. P(x, s¥,)dx,0y, =a (A-4) 


a 


and since negative rounds cannot be fired, 


P(x 9¥,) sane (A-5) 


Equation (A~2) can then be written as 


u(x, s¥,) = ie J P(x, sy,)in q (x,-X,9¥,-¥,) dx, dy, - (A-6) 


The complete problem was minimize 
Oo @® 


~“ul xX 
d= J J PrlX, 94) © mK oVeax, ay, 


subject to 
Co mo 
J :. P(X, 9%) dx, dy, = 1 
oOo 6 
| if es on seb hy = N. 


> = 
P(X V4) 0, Peau) 0 
Weiss continues his solution using techniques similar to those of Morse 
and Kimball described earlier. For any arbitrary P (x, 57.) add a small 


increment 6 to P(x, sy.) over the interval Ax, at x_ and Ay at Yaa ° 


al 
Now examine the change in ¢ The procedure can be illustrated as shown 


in Figure 3. 


ae 











P, (fa,4a) 


Figure 3. Illustration of small increment 6 


added to P(x, »7,) over interval dx at x . and a at Yaa: 


al 


The change in ¢ can be expressed as 


$ —§ p ee = YS 
P(x, »¥,) + AP (x,97,) AM ata : (A-7) 


From equation (A-3), 


o oo : oo , 
= J. J. Pp eV 4) a t l. J. (Py (x17) +4P (x, ¥,) Jl a (4X9 1 4-Vag ) 


: dx_ dy, f 


-exp{- f Jo PyGegsyg)im a(2,-% 1 VE-Vay) dx, ay f wats 


Ae = a J P(X 934) Exp iS - P(x, .y¥,)In q(x,-X5 Y,-Y,7 dx, dy, 


£ fltgeyry i aera reraings ios (LF Nous 
eLn (4a Fe Fa) Va) dx, dy, » 

= J i Pace) | ‘| Py(xa9T, na(x,-% 4 9V4Var aa 
ae dt 1 t _o =~ 


Oo om 


a i AP, (x, 0, )in a (%,-% 7 V4) 9, 
|e ce ces aa | dx, dy. 


2 





CoO © 


co © -u(x, s¥,) “J J OP, (x, .¥,)in a(x,-% 9 s¥4-¥ yy) 
= £ Stare. 
‘dx dy (A-8) 
a~a 
=1 | dx, dy, | 
since P(x, s¥,)=6 and Ln gq (SI VV en) 1s approximately constant over 


a very small area dx dys 


wey -§ px, dy, In a(x,-*,4994-V,7) 


Ad ~ | i) Prtx,»7,Je e -1 |ax,ay,,. (A-9) 


If the bracketed term is now expanded in a Taylor's expansion of e and 
a first order approximation is applied, then 


u(x, 97 


© ) = 
Ag a I P(x, »¥,Je | ~8 dx, by, Ln q(X,-X1 »V4-Va, ) |ox, oy, 


ae u(x, v4) 
~ —6dx Ay, ‘a ie Pr(xpr7,)e Ln a(x,-%, 904 V yg) Wy: 
Weiss defines 
a u(x, 1V,) 
JL Poly ove)e In aya g Vay) Oy Gy4, = pry VQ): 
and then 
Lo. = p(x, 4297) Ax, Ay... (A-10) 


Now either P (x, sy.) 2 On P(x, 5¥,) = 0. That is, either some rounds 
or no rounds are aimed at the immediate vicinity of X97 5° jie 
P(x, sy.) > 0, 5 can be added or subtracted from it; if P(x, ¥,) = 0h. 
§6 can only be added to it. 
: : : S: (ey: 
Weiss considers two points X49% 54 and X609V 50 where P| Conn 0 


6 is added to P(x, sy.) at one point and subtracted from the other. The 


‘ 


change in can be determined from equation (A-10). 
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A? =~A% (x, To4) - A® CT): 
ae es A 
© (x49 44) a © (5 9V go) bax, Ya) 
~| (249% a1)-? (Xgp2%g0) | PES EST (A-11) 


If (A-11) is positive then $ can be reduced further by changing the sign 
of 5, which it is permissible to do since P(x, sy.) > Q at both points. 


It is clear that if 
0 (x0 45¥n4) < —— then A@ <0 ni See 
p (049954) > p (X59V a0) then Aé < 0 ieewOm Oe 
The conditions for an optimum exists when 
0 (x4 5734) = 9S pote = constant C, (A-12) 
for all points where P (x, »¥,) ZeOe 


Weiss then examines a point X39 93 where P(x, 5) = 0. Cae 
can only be added to at this point, Comparing 1 ae to the point 


where DMCS sey = 0, 


A@= E (x, 35Y a3) = o(X, 4 hee) | SAK AY,» 
which is negative if and only if 
< = C. = 
0 (x, 25,3) le) (x, 1¥.4) C (A 13) 
That is to say $ can be decreased only if 0 (X39 a3) <C€. At the 


optimum solution $ can no longer be decreased. 


The original problem is now reduced to solving the following conditions: 


= i iy > 
p (X4,57,4) = C if P,(x,,y,) > 0 fet 
ac if P(x, 575) = 01 


ZS 





ons 


eS = oy) 
ELEM TACs ey) = = J J Pye rVy) BAC, Xa FE Tq IO, 


and u(x, ¥,) is as defined in equation (A-6). 

Weiss has determined necessary and sufficient conditions for the 
aimpoint distribution to be positive. The problem now has been reduced 
to solving equations (A-14) which is a difficult problem and to which 


there exists no general solution in closed form. 
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